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Abstract. We study the interfacial adsorption phenomena of the three-state 
ferromagnetic Potts model on the simple cubic lattice by the Monte Carlo 
method. Finite-size scaling analyses of the net-adsorption yield the evidence 
of the phase transition being of first-order and ksTc/J = 1.8166(2). 
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1 Introduction 



The Potts model is a multi-state statistical mechanical model for spin systems 
0, H|, ||. The Potts spin variable takes a value of 1, 2, ... and g. The Bragg- 
Williams approximation predicts that it exhibits a first-order phase transition 
for q > 2 in all dimensions 0. 

The properties of the two-dimensional models have been well understood. 
According to the duality argument, the transition temperature is known for 
all q exactly [2|]. Baxter || proved the phase transition to be of first-order 
for q > 4 and second-order for q < 4. There is the conjecture for the thermal 
and magnetic exponents y t || and |6|, 0] for q < 4. 

There is no exact result in three dimensions || . There is a problem whether 
the order of the transition of the three-dimensional three-state Potts model is 
of first-order or not. It has been attacked by many authors and the result is 
that the system occurs a weakly first-order phase transition [R, [| [IT], O, O, 



13[ , [LJ, [15|, [16], [17|, |18fl . The transition is characterized by the small jump of 



the energy. 

The specific heat shows the singularity of the delta function type as a 
result of the discontinuity of the energy at first-order phase transitions. The 
behaviour can be seen only in systems with infinite lattice size. In systems 
with finite lattice size the singularity of the energy is smoothed out and that 
of the specific heat is rounded off. The finite-size scaling theory can predict 
behaviour of infinite systems through informations of finite systems |L9| , p0| , |TJ . 

The Monte Carlo method is a useful tool to investigate phase transitions 



and critical phenomena |f22|. p3[ . Monte Carlo simulations are necessarily car- 
ried out on systems with finite lattice size. With help of the finite-size scaling 
theory we can study thermodynamic behaviour of various physical systems. 
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The specific heat of the system with linear size L has a maximum C max (L) 
at a temperature T^^L). The finite-size effects are governed by the dimen- 

5§, |5§: 



sions ti at first-order phase transitions 



C max (L) ~ L° 



T^ ax (L) - T c ~ L 



and by the critical exponents a and z/ at second-order |L9], ^(], 51] 



The difference of the power has been used to determine the nature of transitions 
0- 

As described above the five-state Potts model on the square lattice under- 
goes a first-order phase transition. Peczak and Landau [27| studied the system 
by the Monte Carlo method. They found that it behaved as if the transition 
was of second-order, i.e. the physical quantities showed pseudocritical be- 
haviour, e.g. C max (L) ~ L 109 rather than L 2 . Thus it is difficult to recognize 



the order of the transition. Yamagata and Kasono |28| , however, could found 
the finite-size effects of first-order type in the interfacial adsorption phenomena 



30] of the system by the Monte Carlo method. 
The aim of this paper is to identify the transition of the three-state ferro- 
magnetic Potts model on the simple cubic lattice and to determine the transi- 
tion temperature by analysing the finite-size effects in the interfacial adsorption 
phenomena with the Monte Carlo method. The phenomena in three dimen- 
sions are studied for the first time. In the next section we describe the inter- 
facial adsorption phanomena and then discuss the finite-size scaling theory for 



them in section |3|. We present the detail of our Monte Carlo simulations in 
section [|. In section [| we analyse the Monte Carlo data. A summary is given 
in section ^. 

2 Interfacial adsorption phenomena 



The interfacial adsorption phenomena p9| , p0| have been observed on the multi- 



state models: the Potts models [|^, [33j, the Blume-Capel model |34|, [35 



and the chiral clock model |36|, |37j . We shall explore them for the Potts models. 

Let us consider the g-state ferromagnetic Potts model on the L d ~ l x (L + 2) 
hypercubic lattice. In the <ith direction the fixed boundary conditions are taken 
and periodic boundary conditions are used in the remaining directions. We fix 
the Potts spin variables in two states 1 and m, which takes the value of 1 or 2, 
at the opposite boundaries in the <ith direciton, respectively. The Hamiltonian 
is given by 

* « 

where S is the Kronecker's delta function; J(> 0) is the strength of interactions; 
the summation for is over all nearest neighbour pairs on the lattice; Yl'i 
and J2'l denote the summations over lattice sites which adjoin to the opposite 
boundaries in the cfth direction, respectively. 

In the system Ttn\2) a ^ a very low temperature an interface can appear 
between two phases 1 and 2 (called the (l|2)-interface). It has been found 
that the remaining states 3, 4, ... and q (called the non-boundary states) 
were generated in the form of droplets at the interface near the transition 



temperature Tq p2]| . The interfacial adsorption phenomena are described 
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by the net-adsorption defined BT], |32| as 



W L Ei^EE [{*M) m - (6(<n,n)\ m ] (2) 

i n=3 

where (■ • •) (112) an d {■ • - ) (1 1 1) denote the thermal or Monte Carlo averages for 
the systems Ti.(i\2) and 7~C(x\i), respectively; the summation for i is over all 
lattice sites except those where the fixed Potts spin variables are located. One 
should note that the summation for n is over the non-boundary states. The 
results of Monte Carlo simulations show that Wl(T) has a finite peak near Tq 



on the finite system 31 



3 Finite-size scaling theory 



Yamagata and Kasono |28] have discussed the finite-size effects [32], Q of the 
net-adsorption. Noting the relation between Wl(T) and the interfacial free 
energy, they have shown 

W L (T) « LQ{tL d ) (3) 



and 



W L {T) « L l - d+yh Q(tL 



yt N 



(4) 



at first- and second-order phase transitions, respectively, where Q is a scaling 
function, t = 1 — T/Tq and yt and yh are the thermal and magnetic exponents, 
respectively 0, ||, |6], 0. We omit the correction terms for simplicity here. In 
two dimensions the finite-size scaling (|3[) and (^) of Wl(T) have been confirmed 
by the Monte Carlo method |28|, |32|, One should note in (|3]) that Wl(Tc) is 
proportional to L at first-order phase transitions regardless of the dimensions 
d. 

Yamagata and Kasono [28 have studied the finite-size effects of Wl{Tq) to 
identify the order of transitions. They have found by the Monte Carlo method 



that Wl{Tq) was proportional to L on the two-dimensional five-state Potts 
model. It is clear evidence that the transition of the system is of first-order. 
As described in section [TJ, the finite-size effect of the specific heat of the system 
could not be seen to be of first-order type. It is important to choose physical 
quantities to be analysed. 

We shall study the finite-size effects of the net-adsorption of the three- 
dimensional three-state Potts model on the same lines. We, however, does not 
know the value of the transition temperature of the system exactly ||. We 
note that Wl(T) has a maximum W max (L) at a temperature T m ^ x (L) on the 
system with linear size L. Thus the scaling function Q(x) does so. Let x max 
be a location of the maximum of ft(x) f58|l . x max satisfies a equation 

~ j l^max) 0. 

ax 

We consider the case of first-order phase transitions hereafter. Since it 
is clear that 

_ ~ T^ ax (L) d 

we get a relation 

TZ x (L) = T c -T c x max L- d . (5) 
Since x max is independent of T and L, W max (L) grows with L: 

W max (L) = W L (TZ X (L)) ~ Lt2(x max ). (6) 

We may expect to estimate the transition temperature and to identify the 
nature of the transition by analysing the finite-size effects of the net-adsorption 
obtained by the Monte Carlo simulations. 
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4 Monte Carlo simulations 



We use the Metropolis algorithm |32L ^] to simulate the systems Wnp) and 
W(i|i) given by (|) with g = 3 on L 2 x (L+2) simple cubic lattice (L = 10, 12, 14, 
16, 18, 20, 24 and 30). The net-adsorption is calculated with (§) of d = q = 3. 
We start each simulation from a high temperature with a random configuration 
and then gradually cool the system. The pseudorandom numbers are generated 
by the Tausworthe method |39|. Measurements at a temperature are over 10 5 
Monte Carlo steps per spin (MCS/spin) after discarding 10 4 MCS/spin to 
attain equilibrium. Near T^ X (L) we observe the physical quantities over 10 6 
MCS/spin after 5 x 10 4 MCS/spin for the systems with L = 14, 16, 18, 20, 24 
and 30. We have checked that simulations from the ground state configuration 
and a random one gave consistent results and there was no hysteresis. We use 
the coarse-graining scheme to calculate the statistical errors JI0]. Each run is 
devided into ten blocks and the standard deviations are obtained from the ten 
subaverages. 

As described in the previous section, we want to estimate the maximum 
of the net-adsorption and its position for each lattice size. Since it is difficult 
to get them from raw Monte Carlo data, we decide to adopt a procedure, the 



.B-spline smoothing ||39|| . We fit our Monte Carlo data by the fourth-order 
.B-spline. 

5 Monte Carlo results 

Hereafter, for brevity, the physical quantities are presented in units fee = 1 = 
J. 

Figure 1 shows the temperature dependence of the net-adsorption Wl(T) 
defined by @ with d = q = 3 for various lattice sizes. The net-adsorption has 
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a finite peak. We plot W max (L) versus L in figure 2. It is clear that W max (L) 
is proportional to L. It is consistent with the prediction (^) by the finite-size 
scaling theory. By using the linear- regression from the data W / max (L) with L 
= 18, 20, 24 and 30 we obtain a relation: W max (L) = 0.092(4) + 0.0382(2)L. 
In figure 3, T^ X (L) is shown for the function of L~ 3 . For systems with large 
lattice size it agrees well with (|). We estimate T c to be 1.8166 ± 0.0002 by 
using the linear-regression from the data T^ X (L) with L = 18, 20, 24 and 
30; T^ X (L) = 1.8166(2) - 41(2)L~ 3 . The value should be compared with the 
previous results (see the table |IJ). 

6 Summary 

We studied the interfacial adsorption phenomena of the three-state ferromag- 
netic Potts model on the simple cubic lattice. The net-adsorption was calcu- 
lated by the Monte Carlo method. The finite-size effects were consistent with 
the predictions (f|) and (§) of the finite-size scaling theory for systems with 
L > 18. It is the clear evidence of the phase transition being of first-order. 
The transition temperature was estimated to be 1.8166(2). It is consistent 



with the recent results 16, 18 



The net-adsorption attains to the asymptotic region for L > 18. The 
specific heat, on the other hand, shows the asymptotic behaviour for L > 
30 [|1^]. We succeeded in identifying the phase transition and getting the 
transition temperature from the data with L < 30 since we investigate the 
net-adsorption. Generally we can get good statistics for small systems on the 
Monte Carlo simulations and it takes the long computer time to simulate large 
systems. Thus we need to select the observational quantity carefully. 
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Figure captions 

Figure 1. The temperature dependence of the net-adsorption of the three-state 
Potts model on the simple cubic lattice with L = 10, 12, 14, 16, 18, 20, 24 and 
30. The full curves are obtained by the _B-spline smoothing. As L increases, 
the shape of the curve becomes sharper. 

Figure 2. The size dependence of the maximum of the net-adsorption of the 
three-state Potts model on the simple cubic lattice. The solid line shows 
0.092 + 0.0382L. Errors are less than symbol size. 

Figure 3. The size dependence of the location of the net-adsorption maximum 
of the three-state Potts model on the simple cubic lattice. The solid line shows 
1.8166 -41L" 3 . 
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Table 1: The transition temperature of the three-state Potts model on the 
simple cubic lattice in units k-& = 1 = J. 





T c 


Author(s) 


1 


.8166(2) 


This work 


1 


.816454(32) 


Alves et.al. 1991 


1 


.8164(1) 


Fukugita et.al. 1990 


1 


.8161(1) 


Gavai et.al. 1989 


1 


.81624(6) 


Wilson and Vause 1987 


1 


.81 


Ono and Ito 1982 


1 


.8169(6) 


Knak Jensen and Mouritsen 1979 


1 


.818(3) 


Herrmann 1979 


1 


.818 


Blote and Swendsen 1979 


1 


.827(1) 


Miyashita et.al. 1979 


1 


.787(5) 


Kim and Joseph 1975 
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